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ABSTRACT 

A  certain  conservation  integral,  the  so  called  M-integral, 
recently  exploited  by  L.  B.  Freund  in  the  calculation  of  stress  intensity 
factors  for  cracks  in  homogeneous  elastic  media,  is  applied  to  the 
calculation  of  energy  release  rates  for  interface  edge  cracks.  Specif¬ 
ically  >  for  an  edge  crack  along  the  interface  between  two  elastic 
wedges  of  different  opening  angles  and  dissimilar,  elastic  properties, 
and  that  is  subjected  to  point  loads  at  the  apex,  a  relation  is  derived 
among  the  length  of  the  crack,  the  energy  release  rate  of  the  crack, 
the  applied  loads,  the  wedge  angles  and  the  material  parameters. 
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significance  and  explanation 


The  study  of  the  strength  of  composite  materials  containing  cracks  or 


are  elastic  interface  cracks,  that  is,  cracks  that  have  formed  along  the  inter' 


face  between  two  bonded  dissimilar  elastic  bodies.  An  important  parameter  for 


the  prediction  of  crack  advance  and  ultimate  bond  failure  is  the  energy  release 


rate  of  the  crack,  that  is,  the  rate  of  change  of  the  energy  of  the  body  per 


unit  advance  of  the  crack.  This  report  illustrates  how  a  certain  conservation 


integral  may  be  used  to  calculate  the  energy  release  rate  for  the  2-dimensional 


mixed  boundary  value  problem, 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


ENERGY  RELEASE  RATE  CALCULATIONS  FOR  INTERFACE  EDGE  CRACKS 
BASED  ON  A  CONSERVATION  INTEGRAL 

A.  Nachman  and  J.  R.  Walton^ 

Introduction 

In  a  recent  paper,  L.  B.  Freund  [1]  demon  trated  the  usefulness  of 
the  M-integral  conservation  law  in  the  determination  of  stress  intensity 
factors  for  2-diinensional  cracks  in  homogeneous  elastic  bodies  with  certain 
geometric  properties.  For  a  detailed  discussion  of  the  M-integral  and  its 
special  features  that  permit  the  applications  given  by  Freund,  the  reader 
is  referred  to  that  paper  and  the  references  contained  therein.  Liberal 
use  will  be  made  in  this  paper  of  many  of  the  results  derived  by  Freund. 

What  is  demonstrated  here,  is  that  the  M-integral  can  be  used  to  deter¬ 
mine  the  energy  release  rate  for  certain  interface  cracks  in  much  the  same 
way  as  for  cracks  in  homogeneous  bodies. 

Three  observations  are  needed  for  this  application.  The  first,  due 
to  Smelzer  and  Gurtin  [2],  is  that  the  J-integral  on  a  small  arc  about  an 
Interface  crack  is  equal  to  the  energy  release  rate,  as  in  the  homogeneous 
case.  One  difference  for  cracks  in  homogeneous  bodies  is  that  the  J-integral 
may  also  be  related  to  the  stress  intensity  factor,  while  for  interface  cracks, 
the  J-integral  is  related  to  a  composite  stress  intensity  factor  which 
has  dubious  utility*  (See  Smelzer  and  Gurtin  [2].)  The  second  observation 
involves  the  nature  of  the  far  stress  field  for  bonded  dissimilar  elastic 
wedges,  and  appeals  to  the  analysis  presented  by  Bogy  [3]  for  bonded 
dissimilar  elastic  quarter  planes.  The  third  observation  is  that  the 
integrand  of  the  M-integral  is  continuous  across  bonded  interfaces  lying 
•long  rad:'  .1  lines  of  the  chosen  coordinate  system. 
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Statement  and  Analysis  ot  the  Problem 


We  present  here  an  extension  of  the  result  obtained  by  Freund  for  the 
elastic  half-space  with  an  edge  crack  whose  corners  are  subjected  to  normal 
and  shear  point  loads.  Specifically,  the  M-integral  conservation  law  is 
applied  to  the  2-dimensional  plane  strain  problem  of  two  infinite  isotropic 
elastic  wedges  with  opening  angles  uj'  and  u" ,  respectively,  and  with 
different  elastic  properties  (exhibited  through  E',  v'  and  E",v",  where 
E  and  v  denote  Young's  modulus  and  Poisson's  ratio),  which  are  bonded 
together  along  one  edge  except  for  a  crack  of  length  £.  extending  from 
the  apex,  and  whose  corners  are  subjected  to  a  system  of  normal  and  shear 
point  loads,  given  by  P',  Q' ,  P"  and  Q".  (See  Fig.  1.)  An  application 
of  the  M-integral  conservation  law  yields  a  relation  among  the  parameters 
V,  f.  and  the  energy  release  rate  of  the  crack.  (Symbols 
without  primes  refer  to  both  materials.)  It  is  clear  that  the  com¬ 
plexity  of  the  problem  precludes  the  determination  of  the  energy  release 
rate  by  first  solving  the  corresponding  boundary  value  problem. 

The  M-integral  is  given  by 


M 


J<Hnixi  '  Tkuk,i*i>ds 


(i) 


where  W  is  the  elastic  energy  density,  is  the  unit  normal  to  C 

(which  we  take  to  be  directed  to  the  right  when  C  is  traversed  in  a  given 

direction)  and  is  the  traction  acting  on  the  material  to  the  left  of 

C.  The  stresses  0^  =  °ji  ar®  related  to  the  elastic  energy  W  and  the 

strains  ■  (u^j  +  uj  j)/2  by  2W  ■  °ijUi  j*  where  denote  the 

displacements.  The  traction  is  given  by  -  °ikn1*  The  elastic 

body  is  assumed  to  be  in  equilibrium  without  body  forces,  ie  o  .  ■  0, 

*J  t  J 

and  the  stress-strain  relation  is  given  by  o. .  ■  ^  .  The  conservation 

J  3Gij 

law  for  M  is  that  M  ■  0  whenever  C  is  a  closed  path  surrounding  a 
simply  connected  region  in  the  body. 


Consider  Che  integral 


M  -  H'  +  M" 


-  ](Wnixi  -  Tkuk>lXl)ds 


•  fiW'n'x  -  T'u'  x  )ds  +  [  (W"n"x  -  T"u"  x  )ds 

Jv  i  i  k  k,ii;  -  J  v  ixi  k  k,i  V  ’ 
r*  r" 


where  the  contours  P  and  T"  are  as  indicated  in  Fig.  1  and  T  ■  P  U  P'. 

It  follows  easily  from  the  discussion  by  Freund,  that  there  is  zero  contri¬ 
bution  to  M  from  the  parts  of  T  along  the  crack  faces  and  the  outside 

edges  of  the  wedge.  Moreover,  the  value  of  M  on  the  vanishingly  small 

dP 

arc  around  the  crack  tip  is  4—  ,  where  P(£.)  is  the  potential  energy  of 
the  wedge  with  a  crack  of  length  JL.  (See  Rice  [A]  and  Smelzer  and  Gurtin  [2].) 
Hence,  on  this  small  arc,  M  is  the  product  of  the  crack  length  and  the 
rate  of  decrease  of  the  energy  with  respect  to  crack  length. 

The  contributions  to  M'  and  M"  from  the  small  arcs  around  the 
corner  points  of  the  wedge  follow  from  a  general  result  derived  by  Freund 
for  an  infinite  elastic  wedge  with  corner  loads.  In  particular,  from 
the  small  corner  arc  in  I*'  we  obtain  the  contribution 


(1  -  v'2)  T~  Fa2  Ft2  1 

E'  |_u'  +  8 inti)'  u'  -  sin uP~  _|  * 


and  from  T" 


v2  .  v2  ~i 

E"  J_w"  +  sinuj"  u"  -  sinu"  J 


where,  employing  Freund's  notation. 


F  '  =  -P'cos(w'/2)  -  Q'sin(cj'/2) , 
a 

F  '  *  P'sin(u)'  /2)  -  Q'cos(u)V2), 

F  "  =  -P"cos(uj"/2)  -  Q"cos(u)"/2)  , 

F  "  =  -P"sin (uj"/2)  +  Q"cos(uj"/2)  . 

On  the  bond  line,  x ^  ®  0,  x^  >  £.,  the  contributions  to  M  from  I”  and 
T"  cancel,  since,  among  the  stress  and  displacement  components,  only  0^ 
is  discontinuous  across  the  interface,  and  on  that  line  n^x^  ■  0  (because 
the  interface  is  radial)  and  T^u^  ^  ■  (U1  +  u ^  1^2)xl* 

It  remains  to  determine  the  contribution  from  the  large  arc  for  which 
it  suffices  to  know  the  far  field  solution  for  two  wedges  bonded  together 
with  no  crack  and  with  apex  loads  P'  +  P"  and  O'  -  Q".  From  the  analysis 
presented  by  Bogy  [3]  for  two  bonded  quarter  planes,  it  is  clear  that  the 
far  field  is  radial.  More  precisely,  if  the  elastic  fields  are  represented 

with  respect  to  polar  coordinates  (r,0),  then  roQQ  and  ra  a  vanish 

re 

for  large  r  uniformly  in  8,  whereas  r^rr  does  not.  Consequently,  on 
the  large  arc  we  may  assume  o^g  and  Ogg  are  zero.  It  is  easy  to  see 
that  in  this  case,  to  satisfy  the  equilibrium  equations  and  the  compata- 
billty  equation,  we  must  take  for  the  stress  field 


The  four  constants  A',  B\  A"  and  B"  may  be  calculated  from  the  four 
conditions 


UQ<r,  0  -*■  0-)  -  iig(r,0  Of),  (3) 

u;<r.6  ♦  0  -)  -  u”(r,0  ♦  0+),  (4) 
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f° 

ro'  (r, 
J  rr 

w* 

f° 

J  r°;r(r* 


fOJ 


6)cos(O)d0  + 


ro"  (r ,0)cos(0)d0  -  -(P’  +  P") , 


rr 


0)sin(O)d0  + 


fw"  •• 

L  «rr<r- 


0>sin(0)d0  =  (Q'  -  Q") . 


(5) 


(6) 


-Li 


Equations  (3)  and  (4)  assert  continuity  for  the  displacements  on  the  bond 
line,  while  equations  (5)  and  (6)  express  the  equilibrium  of  tractions  on 
the  wedge  -a)'  <  0  £  w",  0  £  r  £  R,  where  R  is  the  radius  of  the 
large  arc.  A  simple  calculation  shows  that  (5)  and  (6)  reduce  to 

A’ (1-cos  2u)' )  -  B’^u'-sin  2u')  -  A"(l-cos  2w")  -  B"(2a),,-sin  2w")  -  4(Q'  -  Q")  (7) 

-A’(2u>'  +  sin  2ui')  +  B'U-cos  2u)')  -  A"(2u)"+sin  2cjk)  -  B"(l-cos  2w")  -  -4(P'  +  P")  .  (6) 

Substitution  of  (2)  into  the  polar  form  of  the  stress-strain  law 
followed  by  the  application  of  (3)  and  (4)  yields  the  relations 

B'/B"  -  A' /A"  -  (1  +  a) / (1  -  a),  (9) 

•-  re  a  is  one  of  the  two  Dundurs  bi-material  parameters  given  by 

(See  Bogy  [3]  .) 


E'(l  -  y"2)  -  E ”(1  -  v’2) 
E* (1  -  V"2)  +  E"(l  -  V’2) 


for  plane  strain 


E'  -  E" 
E*  +  E” 


for  generalized  plane  stress. 


is  now  an  easy  matter  to  solve  equations  (7),  (8)  and  (9)  for  A’,  B’,  A" 
and  B”.  In  particular,  we  obtain 

A’  *=  (a(Q’  -  Q")  +  b(P'  +  P"))  (4/d) 

B’  =  (-a(P'  +  F")  -  c(Q'  -  Q") ) (4/d) 
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where  a 


A"  and  B"  may  now  be  calculated  from  (9) 


As  observed  by  Freund,  for  the  stress  state  (2),  the  Integrand  of 


Consequently,  the  contribution  to  M  from  the  large  arc  is 


(A'cosd  +  B'sinfl)  a0 


(A"cos0  +  B”sin6)d0 


(A'cos0  +  B'sin0)d0J  (10) 


(b(P'  +  P")  +  c(Q'  -  Q,,)‘t  +  2a(P'  +  P'^CQ’  -  Q"))(2/d) 


Line  (10)  follows  from  (9)  and  the  observation  that 


whereas,  line  (11)  is  derived  by  simple  but 


tedious  algebraic  manipulations 


J 


Combining  the  contributions  to  M  from  the  large  arc,  the  two  small 
arcs  at  the  apex  and  the  arc  around  the  crack  tip  and  appealing  to  the 
conservation  law,  we  obtain 


2 

1  51  =  ~  ~e',V  )  I-CbCP ’  +  p”)2  +  c(Q’  -  Q")2  +  2a(P '  +  P")(q’  -  Q”))(2/d) 


+  ■|(P,2(2ui’  -  sin  2u')  +  Q,2(2id’  +  sin  2io')  +  2P'Q'(cos  2id’  -  l))/(id,2-sin2id’) 


(12) 


+  I(I~  q)(P"2(2u)"-sin  2id")  +  Q"2(2u"  +  sin  2u>")  +  2P"0"(cos  2id"-l)) /(id"2-sin2id")  ] . 


An  important  special  case  of  (12)  is  that  of  two  bonded  dissimilar  quarter 
planes  with  an  edge  interface  crack.  Setting  u'  =  id"  «  tt/2  in  (12)  yields 


dl  V 


UP*  +  P")2  +  (Q*  -  Q")2  -  (P’  +  P”)(Q’  -  Q") )  (1  +  a)/((2o)2-iT2) 


+  2((P'2  +  (fff)P"2)  +  (Q-2  +  (f4-“>Q"2>  -  i  (P'Q'  +  (fr^)P"Q”))/(7r2-4)l. 


(13) 


It  should  be  noted  that  if  P’  -  P"  =  P,  Q'  «  Q"  =  q  and  a  -  0,  then 
(13)  reduces  to  the  result  obtained  by  Freund  for  identical  quarter  planes. 

Of  course,  when  ct  •  0,  the  response  of  the  bonded  quarter  planes  is  the  same 
as  for  a  homogeneous  half-space;  as  described  by  Dundurs  15],  the  two  quarter 
planes  are  "consonant  in  tension  parallel  to  the  interface." 

Another  interesting  case  in  (12)  is  when  a  -  ±1.  Due  to  symmetry  we  consider 
only  a  -  -1  corresponding  to  which  (12)  becomes 


-8- 


i~  -  [P,2(2w'-sin  2a)1)  +  q,2(2u)'+sin  2w')  +  2P,q'(cos  2w,-l))/(u),2-sin2^'  . 


It  should  not  be  surprising  that  in  (14)  only  P'  and  q '  appear,  since 
a  ■  -1  corresponds  to  E"  ■  °°. 

It  should  also  be  noted  that  for  P'  =  P"  E  P,  q'  =  q"  E  q  and 
(j'  =  c  tt/ 2,  the  level  curves  of  l  ^  in  the  (P,q)-plane  are  ellipses 
(straight  lines  if  a  *  0)  centered  at  (0,0).  From  this  we  may  conclude 
as  did  Freund,  that  crack  extension  may  result  when  the  unloading  of  P  and 
q  occurs  along  certain  paths  in  the  (P,q) -plane.  Moreover,  this  obviously 
is  the  case  also  in  (14)  with  u'  V/2.  It  is  evident  in  (12), 
that  in  general,  regardless  of  the  values  of  u'  and  uj"  and  a,  this 
phenomenon  is  to  be  expected. 

The  authors  wish  to  acknowledge  with  much  appreciation  several  helpful 
and  encouraging  discussions  on  this  problem  with  Professor  L.  B.  Freund. 
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